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Abstract 

It is argued that the use of the initial Gaussian energy density profile for 
hydrodynamics leads to much better uniform description of the RHIC 
heavy-ion data than the use of the standard initial condition obtained 
from the Glauber model. With the modified Gaussian initial conditions 
we successfully reproduce the py-spectra, i?2, and the pionic HBT radii 
(including their azimuthal dependence). The emerging consistent pic- 
ture of hadron production hints that a solution of the long standing 
RHIC HBT puzzle has been found. 

1 Introduction 

Relativistic hydrodynamics of the perfect fluid may be considered as the standard framework 
for the description of the intermediate stages of relativistic heavy-ion collisions [1—7]. However, 
despite the clear successes in reproducing the particle transverse-momentum spectra and the 
elliptic flow coefficient V2, the typical approach based on the relativistic hydrodynamics cannot 
reproduce correctly the pion correlation functions. In particular, the ratio of the so called HBT 
radii R out and i? s ide comes out too large, exceeding the experimentally measured value by about 
20-50%. Very recently, we have found [8] that the consistent description of the soft hadronic 
observables may be achieved within the hydrodynamic model if one makes a modification of 
the initial conditions - the initial energy profile obtained in most cases from the optical Glauber 
model should be replaced by the Gaussian profile. 

Our framework consists of the recently developed 2+1 boost-invariant inviscid hydrody- 
namics [9-11] linked to the statistical-hadronization model THERMINATOR [12]. The initial 
eccentricity is obtained from the Monte-Carlo Glauber model GLISSANDO [13]. The simula- 
tions done with GLISSANDO include the eccentricity fluctuations [14-22]. In each simulated 
event the distribution of sources (a mixture of the wounded-nucleon contributions and the binary- 
collision points) is first rotated to the principal-axes frame and then histogrammed. As a result 
one obtains the two-dimensional profile that takes into account the fluctuations of the principal 
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Fig. 1: The transverse-momentum spectra of pions, kaons and protons for the centrality bin c = 0-5% (upper panel), c 
= 20-30% (middle panel), and the elliptic flow coefficient V2 for c = 20-40% (lower panel), plotted as functions of the 
transverse momentum and compared to the RHIC Au+Au data. 

axes with respect to the reaction plane. This procedure determines the initial energy distribu- 
tion in the transverse plane which is parameterized as the two-dimensional Gaussian and used 
as the initial condition for the hydrodynamics. The main characteristic of the hydrodynamic 
stage is the use of the realistic equation of state which interpolates between the lattice simula- 
tions of full QCD and the hadron-gas model. The final stage of the evolution is described with 
the help of the Monte-Carlo thermal model THERMINATOR which simulates hadron emission 
from the freeze-out hypersurface delivered by the hydrodynamic calculation. We assume the sin- 
gle freeze-out scenario [23,24] with the universal final temperature Tj =145 MeV. Besides the 
final temperature our model has essentially two additional parameters: the initial temperature Tj, 
fixing the absolute normalization of the spectra, and the initial time for the start of hydrodynam- 
ics ro =0.25 fm. Of course, for each centrality class we fix the geometric parameters a and b 
(i.e., the widths of the initial Gaussian energy distribution) by the GLISSANDO simulations as 
explained above. 




Fig. 2: The pion HBT radii 7? s idc , -Rout , -Riong, and the ratio -Rout/Rside for central collisions, shown as the 
functions of the average momentum of the pair and compared to the RHIC Au+Au data. 

2 Results 

In Fig. Q] we show our results describing the transverse-momentum spectra of pions, kaons, 
and protons for the centrality classes c = 0-5% and c = 20-30%. Fig. [Qpresents also our results 
describing the elliptic flow coefficients V2 for the centrality class c = 20-40%, plotted as functions 
of the transverse momentum. The spectra and the elliptic flow are compared to the RHIC data [25, 
26]. We observe a very good agreement between the model predictions and the data. The small 
excess of the theoretical proton V2 above the data may be attributed to the lack of rescattering in 
the final state. 

In Fig. [2] we present our results on the pion HBT radii i? s idc . ^out» -Riong, and the ratio 
^out /Aside f° r central collisions, compared to the RHIC data [27]. The left panel shows our best 
results obtained with the traditional Glauber initial condition [28], while the right panel shows the 
results obtained with the Gaussian initial condition [8]. One can see that a very good agreement 
between the data and the theoretical model predictions is achieved in the case where the Gaussian 
initial condition is used. In particular, the ratio R out / R sl d e is well reproduced. We note that the 
calculation of the radii does not introduce any extra parameters. All the characteristics of the 
emitting source were already fixed by the fits to the spectra and V2. 

Finally, in Fig. [3] the results describing the azimuthal dependence of the HBT radii are 
plotted [29]. Here R 2 ((fr) = Rq + 1R\ cos(2</>). Again, we observe a very good agreement 
between the data and our model for different centralities and different average momenta of the 
pion pairs kx- 
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Fig. 3: Results for the RHIC HBT radii and their azimuthal oscillations. For each value of iVp a rt on the horizontal 
axis the experimental points (filled symbols) and the model results (empty symbols) are plotted. The points from top 
to bottom at each plot correspond to kr contained in the bins 0.15-0.25 GeV (circles), 0.25-0.35 GeV (squares), and 
0.35-0.6 GeV(triangles). The top panels show i?o U t,o> -^-side o> an< ^ ^fong o> me bottom panels show the magnitude of 
the allowed oscillations divided conventionally by R^idcfi- 

3 Conclusions 

The results presented above indicate that it is possible to achieve a uniform description of the 
RHIC heavy-ion data collected at the highest beam energies in the soft hadronic sector using 
the hydrodynamics of perfect fluid with the Gaussian initial condition for the energy density. 
In particular, it is possible to describe the transverse-momentum spectra and the elliptic flow 
coefficient t>2 simultaneously with two-particle observables such as the HBT correlation radii 
(the preliminary results show also that the correlations of the non-identical particles are well 
reproduced in our model). Our finding shows that there exists a solution to the long standing 
RHIC HBT puzzle understood as the impossibility of the consistent description of the spectra 
and the HBT radii in a single hydrodynamic approach. 

The use of the Gaussian initial profile leads to a faster development of the initial transverse 
flow which makes the system evolution shorter. At the same time the transverse size of the system 
at freeze-out is slightly larger (as compared to the standard Glauber scenario). These two effects 
put together lead to the desired reduction of the ratio -R ut/-Rsidc- Another important effects 
helping to describe correctly the data consists of the use of the semi-hard equation of state (with 
no soft point leading to the extended duration of hadronization) and of the adoption of the single 
freeze-out scenario which also reduces the emission time, hence, decreasing the radius i? ut- 

Of course, the open question remains to find the microscopic mechanism leading to the 
Gaussian initial conditions. Needless to say, this problem goes far beyond the straightforward 
application of the hydrodynamics that were discussed here. Alternatively, one may think of 
other modifications of the initial conditions, such as introducing the initial transverse flow [9, 



30,31] or separating the system into a thermalized core and an outer mantle/corona consisting of 
independent NN collisions [32-34]. Yet another direction is to study the effects of viscosity [35]. 
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